A RADIATION CONDITION FOR UNIQUENESS IN A WAVE 
PROPAGATION PROBLEM FOR 2-D OPEN WAVEGUIDES 
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Abstract. We study the uniqueness of solutions of Helmholtz equation for a problem that 
concerns wave propagation in waveguides. The classical radiation condition does not apply to our 
problem because the inhomogeneity of the index of refraction extends to infinity in one direction. 
Also, because of the presence of a waveguide, some waves propagate in one direction with different 
propagation constants and without decaying in amplitude. 

Our main result provides an explicit condition for uniqueness which takes into account the phys- 
ically significant components, corresponding to guided and non-guided waves; this condition reduces 
to the classical Sommerfeld-Rellich condition in the relevant cases. 

Finally, we also show that our condition is satisfied by a solution, already present in literature, 
of the problem under consideration. 
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1. The problem of uniqueness for the Helmholtz equation. Let E C 

be a (possibly empty) bounded closed surface. It is well known that the Dirichlet 
problem 



Au + k 2 u = f outside E, 
u = U on E, 



(1.1) 



has not an unique solution. If k = (Poisson's equation), in order to obtain the 
uniqueness, it is required that the solution vanishes at infinity. If k ^ 0, that is not 



sufficient anymore. In fact, there are two different solutions of (1.1) which vanish 



at infinity, representing the outward and inward radiation. Hence, an additional (or 
different) condition at infinity is needed. 
The first condition we can add is that 



lim R 

R. — >oo 



du 

M 



iku 



0, 



(1.2) 



uniformly; this is the so-called Sommerf eld's radiation condition. Here, denotes 



OR 



the radial derivative of u. The physical meaning of this condition is that there are 
no sources of energy at infinity. Moreover, it assures that, far from the surface E, u 
behaves as a wave generated by a point source. 



Stated as in (1.2 1 and together with the assumption that u vanishes at infinity, 



this condition is due to Sommerfeld, see [Sol] and |So2] (see also |Magl| and |Ma g2 ). 
The vanishing assumption on u was dropped by Rellich (see [RelJ, who also proved 



that (1.2) can be replaced by the weaker condition 



lim 

R. — >OG 



dB E 



du 
OR 



— iku 



da = , 



(1.3) 
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where Br is the ball centered at the origin with radius R. In the same paper, Rcllich 
also proved that a radiation condition can also be given in the form 



du 
dR 



iku 



dx < +00. 



(1.4) 



Condition (1.4) can be considered the starting point for our work, as we are going 



to explain shortly. Before describing our results, we cite some generalizations of the 
work of Rcllich. 



When n is a function which is identically 1 outside a compact set, (1.3 1 still 
guarantees the uniqueness of a solution of 



Au + k 2 n(x) 2 u — /, x e 



see [Mil] and [Sc] and references therein. 

Several authors (see |Rel] |Mi2| |RS| |Zh| |PV] jEi]) studied the case in which n 
is not constant at infinity, but has an angular dependency, say n(x) — > n oc (x/\x\) 
as I a; I — ► 00, and it approaches to the limit with a certain behaviour. Among these 
papers, we want to mention the results in [Zh] and |PVj . where the authors proved the 
uniqueness of solutions of the Helmholtz equation by means of the limiting absorbtion 
method and by introducing the radiation condition: 



lim — 

R->+oo R 



du 
dR 



ikn r 



dx = 0. 



Here, the assumptions on n are so that the energy cannot be trapped along any 
direction, but it radiates toward infinity. That is in accordance with |PS] , where 
the authors point out that the Sommerfeld radiation condition, since it involves the 
dimension, is inappropriate for problems admitting a lower dimensional solution (a 
plane wave). 

The present paper is motivated by the study of wave propagation in optical waveg- 
uides. In particular, we shall study the uniqueness of solutions of the two-dimensional 
Helmholtz equation 



where n is of the form 



Au + k 2 n(x) 2 u — /, (x, z) G 



n co (x), \x\ < h, 



\x\ > h; 



(1.5) 



(1.6) 



here bounded function and n c i is a constant; thus, (KL6| models the index of 



refraction of a rectilinear open waveguide of width 2h (subscripts co and cl refer to 
the core and cladding of the waveguide) (see |SLj V 

We observe that functions n of type ( |1.6| are not considered in the works cited 
before. In fact, the most important feature of optical waveguides is the presence of 
waves confined inside the waveguide (guided modes) which are oscillatory and never 
decaying along the axis of propagation (z-axis) . It is easy to show that a pure guided 
mode supported by the Helmholtz equation does not satisfy the radiation conditions 
above retrieved (as already pointed out in |PSj ). Functions n similar to (1.6 1 were 
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considered by Jager and Saito in [JS1 -[JS2 ; however, their assumptions on n do not 
admit the occurrence of guided modes. 

As far as we know, the only works dealing with uniqueness in an optical waveg- 
uide setting have appeared in the Russian literature (see |Rei] |Noj [NS [KNH] and 
references therein). However, the Reichardt condition studied therein only deals with 
guided modes and does not apply to the total field. 

The main result of this paper is Theorem |2.3[ where we present a new radiation 



condition that guarantees the uniqueness of a solution of (1.51 with n given by (1.6 1. 
We observe that, if we suppose that no guided mode is present (this is possible by 
choosing special parameters in the function n), our radiation condition reduces to 



(1.4). In this setting, our results provide a different proof of special cases studied in 



jRelj and |J5T]-|J55|. 

The key ingredients of our proof are essentially four: (i) if (1.5) possesses two 
solutions satisfying our radiation condition, then their difference w must belong to 
the Sobolev space H 2 (R 2 ); (ii) as a consequence of (i), the Fourier transform of w in 
the z-direction (parallel to the fiber's axis) is square integrable for almost all x G R 
and satisfies an ordinary differential equation in x; (iii) the only square integrable 
solution of such an equation is identically zero; (iv) the proof is then completed by 
using an appropriate transform theory in the x-direction and repeating the arguments 
in (ii) and (iii). This scheme will be carried out in [J2] 



In |MS] the authors derived a solutiorj^] for the problem ( 1.5 1-( 1.6 ) in terms of 



a Green's function. Section [3] is devoted to prove that such a solution satisfies our 
radiation condition. This will be done in three steps: in §3.1| we derive a representation 
of the solution as a contour integral; in §3.2| we prove uniform estimates for the non- 



guided part of the spectrum-based solution; in { 3.3 we carry out the proof by testing the 
radiation condition on the guided part and using the asymptotic estimates obtained 
in fHl 

We wish to observe that the results in the present paper can be easily adapted to 
prove the uniqueness of a solution for the Pekeris waveguide problem (see |Wij and 

~])- 

2. A new Rellich-type condition and a uniqueness theorem. In this sec- 



tion we shall state a radiation condition, that generalizes ( 1.4 1, and prove our unique- 
ness result. 



2.1. Preliminaries. We recall the relevant results of |MSj which will be useful 
in the rest of the paper. 



In |MS] a Green's function G for (1.51 is constructed: a solution of (1.5) is 



u{x, z) = J G(x, z; & 0/(f, C)d£dC (2.1) 

R 2 

where 

+ f° J\z-C\ y/k 2 nl-\ 

G(x,z;£,0= E / o- /i 2 2 i v j (x,X)v j ^,X)d Pj (X). (2.2) 

Here 

n* = max n, 



1- We will refer to such a solution as the spectrum-based solution. 
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m=l d2 



VX^d 2 



(A A) 2 + #(M) : 



for all 77 € Cq°(M) (( , ) is the usual dual product), 



+00 



Vj(x,X J m ) 2 dx 



V<P - XL 

h 

Vd 2 -X J m J 4> ] {x,Xln) 2 dx + <pj(h,X J m ) 2 



(2.3) 



Also, Vj(x, X) arc linearly independent solutions of 

v" + [A - q(x)]v = 0, in R, 
where q(x) — k 2 [n 2 — n(x) 2 ], and have the following form: 



(2.4) 



(/)j(h, A) cos Q(x — ft) + ^ 3 q' A ^ sin <5(a; — h), if x > ft, 

if |z| < ft, (2.5) 



Vj(x,X) = { 4>j(x, A), 

4>j(—h, A) cos Q(x + h) + sm q( x + ft), if x < —ft, 



for j = s,a, with Q = VA — d 2 , d 2 = k 2 (n 2 — n 2 cl )\ the (f>j's are solutions of (2.4 1 in 
the interval (— ft, ft) and satisfy the initial conditions: 

0,(O,A) = 1, ^(0,A) = 0, 
6,(0, A) = 0, #,(0, A) = n/A. 

(The indices j = s,a correspond to symmetric and antisymmetric solutions, respec- 
tively.) 

We notice that (2.2 1 can be split up into two summands, 

G = G 9 + G rad , 

where 

M j i\z — £\y/k 2 n 2 — \ j 

G 9 (x,z;£,C) = E E 6 / *= Vj(x^in>j(t^l)rL (2-6) 

je{s,a} m=i 2iy fc 2 n 2 — Am 



and 



with 



+00 



j£{s,a} d 2 



x'^- x ^-\x >■ (2 - 7) 



X-d 2 



(x-d^^xy + Mh.xy 



(2.8) 



G 9 represents the guided part of the Green's function, which involves the guided 
modes, i.e. the modes propagating mainly inside the waveguide; G rad is the part of 
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the Green's function corresponding to the non-guided energy, i.e. the energy radiated 
outside or vanishing along the waveguide, which we denote by 

u rad = J G rad (x,z;Z,C)f(t;,0- (2.9) 

K 2 

It exists a finite number of guided modes, which corresponds to the finite number 
of roots of the equations 

VtP-A^tMj+^CM) =o, je{s,a}, 

laying in the interval (0, cP). We shall denote by X J m , m — 1, . . . , Mj, j — s, a, these 
roots. Each Vj(x, X J m ) decays exponentially for \x\ > h as it is clear from the formula: 

{^■(/i,A4)e-Vd 2 -^(^-ft) ; x> h, 
^'(^A4), \x\<h, (2.10) 

ct> J {-h 1 \i n )e^ d2 - xL{ - x+h \ x<-h. 



We notice that G 9 is bounded and oscillatory in the z direction, because \j k 2 n 2 — XI 
is real for every m = 1, . . . , Mj, j = s, a. 



Remark 2.1. The functions o~j(X), j £ {s,a}, given by (2.81, are meromorphic 
functions of X S C, real-valued for A 6 R and with poles that are real and simple (see 
\CL1 .\Ti j ), which corresponds to the values X 3 m , m = 1, . . . , Mj, j — s, a. 

To simplify notations, we shall denote by 7;, I — 1, . . . , M, M = M s + M a , the 
values X 3 m , m = 1, . . . , Mj, j — s, a, ordered according to the natural ordering on the 
real line, and by 7* their maximum. With these premises, we set 

e(z,70=^4- ^ 



From (2.6 1 and (2.3), it is clear that the guided part G 9 can be written as 

M 

G 9 (x,z;t,() = Y,G 9 i(x,z;t,(), 
1=1 

where 

G 3 ,(x, z; £, C) - £ -^e(x, 7l )e^, 7l ), (2.12) 

with 

A = Vk 2 n 2 - 7J . (2.13) 

Let s £ R; we will denote by L 2 ' S (K 2 ) the weighted Lebesgue space consisting of 
all the complex-valued measurable functions u such that (1 + x 2 + z 2 ) s \u{x, z)\ 2 is 
summable in M 2 , equipped with the natural norm defined by 

\ u \i s = I \u(x,z)\ 2 (\ + x 2 + z 2 ) s dxdz; 
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L 2 ' S (M 2 ) is commonly used when dealing with solutions of Hclmholtz equation (see 
|Ag| and |Le] ) . In CM we proved that the spectrum-based solution (2.1 1- (2. 2 1 derived 
in |MS] belongs to £ 2 ' S (IR 2 ), for s < -1, if / £ L 2 <~ S (]R 2 ). 



The following lemma will be useful in the next subsection. 
Lemma 2.2. For s < — 1, let w e L 2 ' S (R 2 ) satisfy 

Au + k 2 n(x) 2 u = 



m 



where n is given by (1.6). Then 

Km u{x,z)e-\ x ^ d ' 2 -'i* = lira u x (x, z) e -^\ /d2 -"" = 0, 

|x| — »+oo |a;| — >+oo 



where 7* = max 7;. 

1</<M 



(2.14) 



(2.15) 



Proof. Since u is a solution of (2.14), from Lemmas A.l and A. 3 in |CM) . we infer 
that both (1 + x 2 + z 2 ) s \ \7u(x, z)\ 2 and (1 + x 2 + z 2 ) s \V 2 u(x, z)\ 2 are summable in 
K 2 . Thus, it easily follows that the function 

^(x, z) = (1 + x 2 + z 2 ) s/2 u(x 7 z) 

belongs to the Sobolev space W 2,2 (M. 2 ). The Sobolev Imbedding Theorem (see Theo- 
rem 4.12 in |AFj ) implies that ^5> £ L°°(IR 2 ) and hence the first limit in (2.151 follows 
at once. 

A straightforward computation shows that \& satisfies the following equation 

A* + b ■ V* + c * = 

in M 2 , where 



b(x, z) 



2s(x, z) 
1 + x 2 + z- 



:, c(x, z) = k 2 n(x) 2 — s 



2 - s{x 2 + z 2 ) 
(1 + x 2 + z 2 ) 2 ' 



Since * £ W 2 ' 2 (M. 2 ), by Theorem 8.10 in [GT], we have that * £ W 3 ' 2 (H + ) where 
H + = {(x,z) £ M 2 : x > h}. Again, by the S obole v Imbedding Theorem, |W| is 
bounded in H + and hence the second limit in (2.151 holds as x — * +00. The same 
limit as x — * —00 holds by a similar argument. □ 



2.2. The radiation condition and uniqueness theorem. We consider a so- 

,...,M, (2.16) 



lution u of (1.5) and define 

ut(x,z) = e(x,ji)U(z,ji), I 



with e(x, 7z) given by (2.111 and where 

00 

U(z, ll )= J M (e,z)e(£, 7 zR, l=l,« 

— 00 

The remainder part of u is 

M 

uq(x, z) = u(x, z) — ui(x, z). 



M. 



(2.17) 



(2.18) 



1=1 
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We introduce a one-parameter family of sets flu, R > 0, such that dfla are level 
sets of a convex and coercive function d(x, z), i.e. f2# = {(x, z) G M 2 : d(x, 2) < i?}. 
With these notations, we state our radiation condition for a solution u of ( |1.5| : 



M 

E 

;=o 



9i> 



(tfeLR < +00, 



(2.19) 



OQr 



with /3 = kn c i and /3;,Z = 1, . . . , M, given by (2.131. 

Notice that, when n = 1, we can choose flu = Br and (2.191 reduces to (1.4), 
since in such a case the guided components are not present. 

The main result of this paper follows. 



Theorem 2.3. There is at most one solution of (1.51 that satisfies (2.191 and 
belongs to u e L 2 > S (M 2 ), s < -1. 

Remark 2.4. vis if will be clear, it is not necessary to specify further the sets 



$Ir in (2.191 to get uniqueness of a solution of (1.5 1 . That means that Theorem 2.3 



holds for any choice of one-parameter family of sets ilji satisfying the above mentioned 
assumptions. 

Of course, a solution of (1.5| satisfying (2.191 may not exist for an arbitrary 
choice of the sets £Ir. In §3, we shall choose a special family of sets £Ir and prove 
that, with this choice, the solution of (1.51 given by ( 2. 1 ) -( 2.2 1 satisfies (2.191. 

We also notice that it is not necessary to choose the same sets in each ad- 
dendum in (2.19); Theorem 2.3 still holds if we replace (2.19) by the more general 

dui 



radiation condition 



M 

E 

1=0 



dv 



ifiiui 



didR < +00, 



(2.20) 



where fi£\ 1 = 0,1, 



,M, are one-parameter families satisfying the above mentioned 



assumptions. 

Theorem 12.31 is based on Lemma \2. 51 and Theorem 12.61 below. 



Lemma 2.5. Let j3 £ K and u be a weak solution of (2.141 Then 

d£, 





dw 




dw 




1 


Q IPW 

av 


W( 


dv 


2 + /? 2 M 2 ) 



(2.21) 



on 



on 



for every fl C K 2 bounded and sufficiently smooth. 

Proof. Since u is a weak solution of (2.14), by Theorem 8.8 in |GT] , we obtain 
the necessary regularity to infer that 



on a 

|2 ;„2 /„, n2i„ .12 



{IV 



u 



div(uVw) dxdz 
k 2 n(x) 2 \u\ 2 } dxdz. 



J {|Vu| 2 + uAu} dxdz 



n 
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Therefore 



Im f w ^ dl = 0, 
J ou 

an 



which easily implies (2.211. □ 

Theorem 2.6. Let u e L 2 ' S (]R 2 ) be a weak solution of (|2.14[) satisfying (|2.19 1 



Then 



M +°° 
^ l dR l 

l -° o aojr 



9i/ 



A 2 M 2 



< +oo, 



and, m particular, 



\uA dxdz < +oo, 



(2.22) 



(2.23) 



£ 2 



/or every 1 = 0,1,..., M . 



Proof. By Lemma 2.5 



Suppose / > 1. Since 



it is enough to prove that each ui, I = 0, 1, . . . , M, satisfies 



(2.141. Then, (2.221 and (2.231 will follow from (2.211 and (2.191 



Vm • dxdz + k / n(x) utp dxdz = 



(2.24) 



for every ip G Hq(M. 2 ), we choose y>(x, z) = e(x,^i)r)(z) with 77 S Cg(IR), and obtain: 

)e'(x, ji)rj(z) + u z (x, z)e(x, ~fi)rf{z)] dxdz+ 



k / n(x) u(x, z)e(x, r yi)n(z)dxdz = 0; 



an integration by parts and Lemma 2.2 then give 

u(x, z)e" (x,ji)rj(z)dxdz — / u z (x, z)e(x,'yi)ri'(z)dxdz+ 



k 2 / n(x) 2 u(x, z)e(x, r yi)n(z)dxdz = 0. 



Since e(x,ji) satisfies (2.4 1, we obtain 



u z (x, z)e(x,ji)n l (z)dxdz + (k n 2 , — ji) / u(x, z)e(x, ji)i](z)dxdz = 



and thus, from (2.171 



U z (z,-fi)r]'(z)dz+ (k 2 n 2 - 7J ) / [/(z, 7i )?Kz)cfe = 



UNIQUENESS FOR 2-D WAVEGUIDES 



9 



for every rj e Cq ( 



Together with (2.4). this formula implies that each ui(x, z), I = 1, . . . , M, given by 



(2.16 1, is a weak solution of (2.14 1. In fact, for tp(x, z) — ip(x)r](z) with i/),r) € Cq ( 



integration by parts gives 



Vu;(x, z) • V<p(x, z) dxdz + k 2 I n(x) 2 ui(x, z)ip(x, z) dxdz 



{e"0,7/) + [7z - l{x)]e(x,^i)}ip(x) dx\ U(z,ji)i](z)dz + 



the same conclusion holds for any y> € Cq(M 2 ), by a density argument. 

Since u and u/, Z = 1, . . . , M, now satisfy (2.14), the same holds for u . Thus, as 



already mentioned, we can apply Lemma |2.5| to each ui, I — 0,1, ... , M, and obtain 



M 

E 



dui 
dv 



an r 



M 

E 



for every i? > 0, and then, since u satisfies (2.191, we get (2.22) and (2.231. □ 



2.3. Proof of Theorem 12. 31 Let u\ and U2 be two solutions; u = u\ — U2 



satisfies (2.141 and (2.191 



From Theorem 



2.6 



we have that u E L 2 (M. 2 ) and, by using Lemmas A.l and A. 3 
in |CMj . we get u e H 2 (R 2 ). Therefore, u(x, •) € L 2 (R) for almost every 16R, and 
the same holds for u x {x, ■) and u xx (x, ■). Hence, we can transform (2.141 by using the 



Fourier transform in the z-coordinatc, 



+00 



u(x,t) = / u(x,z)e lz dz, for a.e. i£l, 



and obtain: 

u X x( x , t) + [k 2 n(x) 2 — t 2 ]ii(x, t) — 0, a.e. x £ 
From Fubini-Tonclli's theorem, the integrals 



-foo +00 +00 +00 

\u(x, t) \ 2 dxdt, dt \u(x,t)\ 2 dx and dx \u(x,t)\ 2 dt 



(2.25) 



-00 —00 



-00 —00 



have the same value, finite or infinite. 

Since u(x, ■) belongs to L 2 (M.) for almost every the same holds for u(x, •) 

and, furthermore, we have 



+00 +00 

|2jj. o_ / l„.f™ ~\|2, 



|u(x,£)| dt = 2n / \u(x, z)\ dz a.e. x£ 
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By integrating the above equation and using Fubini-Tonclli's theorem, we obtain 

+00 +00 +00 +00 

J \ii(x,t)\ 2 dxdt — J dx J \u(x,t)\ 2 dt = 2tt J dx J \u(x, z) \ 2 dz 

R2 —00 —00 

= 2tt \u(x, z)\ 2 dxdz < +00. 



-00 —00 



Therefore u(-,t) £ L 2 (R) for almost every tel. 
From ( |2.25[ ), it follows that 

u(x, t) — a(t) cos \J A — d 2 (x — h) + b(t) sin yA — d 2 (x — h), for x > h, 
where A = k 2 n 2 — t 2 and d 2 = k 2 (n 2 — n 2 t ). Since 



+00 



\u(x, t) \ 2 dx > / \u(x,t)\ 2 dx, 



we obtain that u(x, t) can be not identically zero only for some values t = A^ G (0, d 2 ] 
and, furthermore, in that case 

u{x,t) = a{t)v s {x,\™) +b(t)v a (x,\™). 

Hence, for some A, B £ R we should have 

u(x, z) = AZ s (z)v s (x, XT) + BZ a (z)v a (x, A™), 

where Zj(z) — e ±z V fc2 ™?- A ™ ; because u is a solution of (2.141. Since u(x, •) € L 2 (M), 
then both A and B must be zero and hence u = on IR 2 \ 

3. The spectrum-based solution satisfies the radiation condition. It will 
be useful to introduce the following function 



x + h, x < — h, 
[•'•J/, = <0, -h<x<h, 
x — h. x > h. 

This section is devoted to the proof of the following result. 



(3.1) 



Theorem 3.1. Let f £ L 2 (IR 2 ) be such that / = a.e. outside a compact 
subset ofM. 2 . Then, the spectrum-based solution (2.1) of (1.5 I is the only solution in 
L 2 ' S (R 2 ), s < -1, such that 



du 2 

— lP U C 

av 



on R 
where O is given by 



M 

1 = 1 n 



dui 
dv 



d£dR< +00, (3.2) 



0Q R 



Q r = {{x,z) £ R 2 : [x] 2 h + z 2 < R 2 } 



(3.3) 
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Figure 3.1. The set U R . Figure 3.2. The contour C. 



(see Fig. 3.1) and Qr = {(x,z)em: 2 \x\,\z\<r} 



Remark 3.2. At the cost of extra computations, it may be proved that Theorem 
3.1 also holds if we replace (2.201 by the more compact condition (2.19) with given 



by ^3|. 

We shall break the proof of Theorem |3.1| up into three steps. First, in |3.1| we 
will derive a handier representation of the radiating part G rad of the Green's function, 
as a suitable contour integral (see Lemma 3.3 1. Then, in i 3.2 we will prove a uniform 

asymptotic expansion for the quantity — if3$G rad on the sets dfln. Such an 



dv 



expansion will be used in ^ 3.3 to carry out the proof of Theorem 3.1 where we also 
test the radiation condition on the guided components of u. 

3.1. Representing G rad as a contour integral. We introduce the following 
functions: 



x - [X\h, 



with [x]h given by (3.1 1, and, for r G 



0;.({xH,d 2 + T 2 ) 



$j(x,T) = (/> 3 ({x} h ,d 2 + r 2 ) + \ j e { Si a}. 



(3.4) 



With these notations, ( |2.5[ ) and (2.8 1 take the more compact forms: 
V:j (x,d 2 + T 2 ) = U$j(x,T)e iT W h +Q j (x,-T)e- iT[x]h } 



and 



for j € {s, a}. 



a 3 (d 2 + T 2 ) = 



1 



$y(/l,T)$j(/l,-T) 



(3.5) 



(3.6) 



Lemma 3.3. Let C be the contour from — ^ + i ■ oo to ^ —% ■ oo shown in Fig. 3.2 
and let G rad be the function in (2.7 1. Then, 



je{s,a} c 



A+(x,Z;t)e 0oa + {x ' z *^ + A-(x,Z;t)e l P oa -( x ' z *>^ 



dt 
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with 



and 



± 1 $j(s^ sin (g,±/? sin t) 

j[X ^ 1 ' 8tH ' $j (h, fa sin t)®j(h, -0 O sin t) ' 



a±(»)2;C,C;t) = ± [£]/,) sin i + |z - C|cosf, 



i G C, and where <frj, j £ {s,a}, is given by (3.4 1. In particular, the following 
equivalent expression for G rad will also be useful: 



(3.7) 



where 



G rad = J gfafrt) e ^o([x]hSint+|z-C|coBt) dt) 

c 

g(x,£,t)= £ [A+(x.,^t)e^ sint + Aj(x,^t)e 

je{s,a} 



i[?]h sin* 



(3.8) 



(Notice that g does not depend on x for \x\ > h.) 

Proof. We first take ( 2.7 1 and make the change of variable r = ^/X — d 2 to obtain: 



G 



rad 



47TJ 



-«j(a;,T 2 + d 2 )vj(£,T 2 + d 2 )o-j(T 2 + d 2 )dr; 



j'G{s,a}_ 

here, we also used the fact that all the relevant quantities subject to integration are 



even functions of r. With the help of (3.5) and (3.6 1, and simple manipulations, we 
can infer that 



G 



rad 



-. £ 



+00 



8ni 



j<£{s,a} 



$j(a,T)$,-(g,T) i[T([x] h +[£] fe )+| z -ClV^F^] 
rJ^-^-r) 

^i(x,r)^(e,-r) i[ T ([,] h -[a h )+|,-ClV^^] 

The conclusion is then readily obtained by splitting the interval of integration up into 
the three intervals (— oo, — /3 ), [— /3o>A)] and (/3q,+oo) and by subsequently making 
the change of variable r = (3 Q sin t, with t EC. D 



Lemma 3.4. For every £, £ fixed, we have: 

QQrad 



dx 



= i/3 / ff (/isigna:,C;t)sinte l/3o([:ElhSillt+|z " c|cos * ) di, 



(3.9a) 



for \x\ > ft. and z ^ £; 



0o sign (* - C) / gix^i^coste^^^+^-^^dt, 



(3.9b) 
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for z^(. 



In particular, on the set (0, C) + given by (3.3 1 7 we have: 

QQrad 



dp 



if3 G rad = O / g(x,£;t)[cost- l}e ll3oRcost dt, 



(3.10a) 



for z ~ Q — R and \x\ < h, and 

QQrad 



dv 



i(3 G rad = i(3 / g(h,^t)[cos(t - -d) - l}e^ oRcos ^-^dt, (3.10b) 



where v is the normal to (0, C)+dfln and we have set [x]h = Rsiwd and z~Q — i?cosz? 
with § S [0, tt/2) and R > h. 

Formulas analogous to (3.101 hold for the remaining values of § in [— 7r,7r). 

Proof. Since z ^ ( and Im([a;]^ sint + \z — £| cost) — > +oo as t — > oo on C, the 
integrands in ( 3.9a I and (3.9b I vanish exponentially as t — > oo on C, since g is bounded 
(see Lemma A.3l. Thus, (3.9a I and (3.9b I follow from an application of Lebesgue's 
dominated convergence Theorem. □ 

QQrad 

3.2. Uniform asymptotic estimates for — i(3 G rad . Aiming to esti- 

a ov 

mate, as R — > oo, the function i/3 G rad given by (3.10), we need to deform 

the contour C to a more convenient one. 

Without loss of generality we can assume that i? € [0, tt/2]. We define the new 
contour C$ (see Fig. 3.3 1 as follows: 



C,9 - |J Tj, 



where 



Si = arccos 



2/?o 



y/WS + <P ~ lAt'' 

(notice that cos Si cosh S2 = 1) and 



$2 = arcsinh 



1M 



2/3o 



ri 


= {t = h 


TT 

+- H2 G C : Rc(cost) = 1, Im(cost) > 0, — — 


<h<-Si, t2>5 2 }, 


r 2 


= {t e C : 


-Si <ti< -Si + <&,t 2 = S 2 }, 




r 3 


= {teC: 


Re[cos (t - 1?)] = 1, Im[cos (t - 1?)] > 0, \ti 


-$\<Si, \t 2 \ <<5 2 }, 


r 4 


= {teC: 


Si + 1? < h < it- Si, t 2 = -S 2 }, 




r 5 


= {teC: 


7T 

Re(cost) = -1, Im(cosi) > 0, - < ti < ir - 


-Si,t 2 < -S 2 }. 



This choice of C$ is suggested by the following three remarks: 

(i) C U Cg does not contain in its interior the poles of g (which correspond to the 
guided part (2.6 1 of G); 

(ii) r 3 is part of the steepest descent path of cos(i — t9); 



(iii) ri,r2,r4,r5 are chosen to complete the contour CUC§ and to fulfill Lemma 3.5 
below. 
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Figure 3.3. The contour C. 



By (i), it is clear that we can write 

QQrad 



. hi -t(3 G rad = z/3 / g(x,Z;t)[cost- l]e^ Rcost dt, 



Co 

for \x\ < h, and 

QQrad 



if3 G rad = ip / g{h, f; i)[cos (t - ■&) - 1] e tl3oRcos ^~^dt, 



dv 
for x > h. 

Lemma 3.5. Let (x, z) £ (0, £) + 8CIr. The following estimates hold for R — > oo: 



dG 



rad 



ip G rad = ip J g(x,£;t)[cost- 1} e^° Rcost dt + 0(e- c(3 ° R ), (3.11a) 



dv 

r 3 



for \x\ < h, and 

QQrad 



dv 

r 3 



ip Q G rad = ip / g{hX]t)[cos{t - d) - l]e lfioRcos(t -^ dt + Oie'^ 011 ), 

(3.11b) 



/or x > h, -d e [0, 7r/2], w/iere 



d 2 -^M . (, \fd 
„ = • mm 1 , 



2 - 7m 



4/3 2 + d 2 - 7M I ' 2P 
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Proof. We shall prove only ( 3.1 lb I since ( 3. 11a I follows analogously. We prelim- 
inarily observe that 



Imcos(£ — #) > c, 



(3.12) 



for t £ ri,r 2 ,r 4 ,r 5 and i? £ [0,7r/2]. From Q3.12|), we easily obtain that 
g(h,£;t)[cos(t-tf) - l] e ^« cos 



Ktt 



< ^(cosh5 2 + l)e- c ^, j — 2, 4, 



where if is a bound for g (see Lemma A.3l. Thus, it remains to prove that 
J g(h,Z;t)[cos(t-$)-l}e lf3oRcos ( t -^dt = 0(e~ cf)oR ), j = l,5, 

uniformly as R — > oo, for $ € [0,7r/2]. We carry out the details for j — 1, the case 
j = 5 is completely analogous. We first use Lemma A. 2 to write that 

£; t)[cos (t - i?) - 1] e *Po*«»(*-*)dt = J(R) + 0{e- c/3oR ), 



since (3.121 holds; here, 



J(i?) = »jSb / 1 + 



2(3q sini 



p(y)d» J [cos (t-tf) - i] e m[R™s(t-i>)+(.h-$) s int] dt> 



with p(y) = d 2 — q(y). Let ip(t) = Rcos(t — i?) + (h — £) sini and 6 = 6x + iS 2 ; an 
integration by parts yields 



V>'((5) V 2/3 sin<5 



p(y)rfy [1 - cos (5 - i?)] 



i(t-t?)-l] 



1? 



#(*) 
2/3 sini 



P(y)dy 



+ (/i- - sini?) 



iip'(t)(cost — cos$) 
2/? sin 2 £ 



From (3.12) and since 



sinh£2 < | cos 1 1 , | sin £| < cosh £2, 
1 



< A)(# + |/i-£|)cosht 2 , \ip'(t)\ > -p Rsinht 2 , for > 2\h - £| coth5 2 , 



for f S Ti, we obtain that J(i?) = O(e- c/3oit ), as R -> 00. □ 
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Theorem 3.6. On OCIr, we have 

QQrad 



dv 



1P0G 



rad 



o 



(3.13) 



uniformly as R — > oo . 

Proof. First, we estimate the left-hand side of (3.131 on the sets (0, C) + diln. 
By Lemma 3.5 we only need to estimate the first addendum in (3.11). We prove the 
estimate for (3.11b I; the estimate for (3. 11a I follows exactly in the same way. 

Since F 3 is part of the steepest descent path, the steepest descent method (see 
|BOJ) suggests to change the variables in the first addendum in (3.11b): by setting 
cos(i — ■&) — ! + iy 2 , we obtain 



-l(h,e,t){co S (t-$)-l}e l P° Rms ( t -^dt = -4ie^ R / yV^ 2 ^MMldy, 



with 2/0 = (sin 5i sinh 62) 2 ■ Thanks to Lemma 



A.3 



g(h, £ t)[cos (t - e *fhR<*»(t-*)dt 



(la 



< 2V2K I y 2 e-^v dy < Kj^ R-i 



where K is a bound of g. Therefore, (3.11b I implies that 

QQrad 



i(3 G rad 



r 7T „ 3 



on the sets (0, () + dflR- 

By using exactly the same argument as before, we can prove that the derivatives 
of G rad are Q(i? ") on the sets (0,0 + dtt R , uniformly as R — * 00; we reach the 
conclusion (3.131 by observing that Vdn R ~ ^(o,C)+an B = (^(R^ 1 ), as i? — > 00. □ 



3.3. Proof of Theorem [37T} Since / G L 2 (E 2 " 
from Corollary 5.1 in [CM] we have that u G L 2 ^ S (IR 2 ), 



and / has compact support, 
s < -1. 



Thus, it remains to prove that (2.1) satisfies (2.191. In order to do it, we shall 
check the following facts: 



(i) if u is given by (2.1 ) and Ui, I = 1, . . . , M, is computed via (2.16 1, the remainder 
part u of u, given by (2.18), equals the function u rad in (2.9|; 



(ii) u satisfies (3.2 1. 

We preliminarily notice that 



+00 



dQr 



du 
dv 



i/3 u 



dldR < +00 



is easily verified thanks to Theorem |3.6| 

The following property of orthogonality is useful to check (i) . 
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Lemma 3.7. Let e(x, 7;), 1 = 1,... ,M, and Vj(x,X), j € {s,a}, be the solutioi 



of ( 2.4 1 given by ( 2. 11[ > and (2.5 1 7 respectively, with A > 0. If A 7^ 7;, £/ien 



+00 



e(x, ji)vj(x, X)dx = 0. 



Proof. We multiply the following equations 

e"0, 7;) + hi ~ q(x)]e(x, 7/) = 0, 
v (x, A) + [A — q(x)]v(x, A) = 0, 

by v(x, A) and e(x,7;), respectively, and integrate in x over an interval (a,b). An 
integration by parts gives: 

6 6 
(7/ -A) / e(x,-yi)v(x, \)dx = / [e(x, ji)v"(x, A) - e"(x, Ji)v(x, A)] dx 



= [e(x, 7z)v'(x, A) - e'(x, 7;)w(x, A)] a . 

The conclusion follows by observing that e(x,7;) and its first derivative vanish expo- 
nentially as |x| — -> 00, while v(x, A) and v' '(x, A) are bounded. □ 



Now, by (2.1 1, (2.161 and Lemma 3.7 we have that 



ui(x,z)= Gf(x,z;Z,0f(Z,0d(id<;, 1 = 1 M, 



with G? given by (2.121 and thus uq = u 



,rad 



To complete the proof it remains to check (ii) for I = 1, . . . , M. When z is large 
enough, we have 



dui 
dv 



ipm = o, l = i,...,M, 



on dQn n {(x, z) : \z\ = R}, since ^ = ±gj- Thanks to (2.10), we easily find that 



dui 
dv 



= 



,-Ry/d 2 —n 



as R — ► 00 on 8Qr n {(x, z) : \x\ = R} and thus we obtain that 

+00 



dui 
dv 



ifiiui 



dQ R 

which completes the proof. 



dldR < +00, l=l,...,M, 



Appendix A. Asymptotic Lemmas. In what follows 
space of functions with bounded variation. 

Lemma A.l. Let T be a non-negative number, q € BV(M) and 

p(x) = d 2 - q(x), i£i 



denotes the 
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Then, the following asymptotic estimates for the functions <5> s and <I> a given by (3.4 1 
hold uniformly as \t\ — > +oo, for i£R and | Imr| < T: 



<f>a(x,T) 



{x 






"hi 


p(y)dy 




L 










Vt 2 + d 2 




I 


IT 






(A.l) 



jr {x}h+a (^_y (A-2) 



Proof, (i) First, we prove an estimate for <pj, j £ {s,a}. From (2.4|, we know 
0;'(y,A) + [r 2 +p(y)]^(y,A) = O, J6[-M], 



that 4>j satisfies 



We multiply the above equation by sinr(x — y), integrate by parts twice and obtain 
the following integral equation: 



<t>j(x,X) 



^(0,A) . 



1 



sin tx + <j)j(0, A) cos tx / p(y)sm.(r(x — y))(fij(y,X)dy. (A. 3) 



We set j]j(x,X) — sup \4>j(s, A)|. Since | sinrxj, | cosra;| < coshr 2 a; (r 2 = Imr), 

s£[0,x] 

from the above equation we have that 



r}j(x,X) < 



\m A)| 



|^(0, A) | 



coshr 2 x + / P(y) coshT 2 (x - y) r)j(y, A) dy, 



and, by Gronwall's Lemma (see |SCj ), we get 
1^(0,A)| 



r)j(x,X) < 



|^(0, A) | 



e ° x 



- f p(y) cosh r 2 (x-y) dy 

x < 1 + r 2 / e sinh r 2 s c?s 



Since < p(y) < d 2 , we have that 



J7j(i,A) < 



1^(0, A)| 



|^(0, A) | 



coshr 2 a; exp 



If we assume |r| > d and x £ [—h, h), we finally get 

d sinh T/i 



A) | < v^cosh r 2 h 



exp 



T 



d 2 sinh r 2 a 
|r|r 2 



j G i>,a}- 



(A.4) 



manipulations we have: 



(ii) Now we prove ( A.l ) and (A. 2 1. Let q £ C 1 (IR). From (3.4 1, by straightforward 
8: 

&Ax,t) -ir^j{x,T) = -p(x) <t>j{x,X); (A.5) 

T 
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by multiplying the above equation by e lTX , integrating by parts twice and observing 
that 



2 / e- i Typ(y)4> j (y,X)dy= / e^p^^y, X)dy 



+ - f e-^yp , (y)cj )j (y,X)dy^ 

IT J 




-iTy 

- — p{y)<t>j(y,x) 

IT „ 
y=0 



it follows that <f>o satisfies 



^{x,r)e- iTX = %(0,r) + - / p{y)^{y , X)dy 



2r 2 



e- iry p'(y)cj>j(y,\)dy-e- i ™p(x)<l) j (x,\)+p(0)<t>^,X) \ • (A.6) 



By setting Mj(x,r) = sup \$>j(s, A)e lrs | and from (A.4|, we get 

se[o,x] 



M,-(x,r) < |^(0, r)| + ^ / p(y)M j (y,T)dy 



+ ^\ C J e T2y \p'(y)\dy + Cd 2 e^ +p(0)\M0^)\ \ , 



for |t| > <i, where C is the right-hand side of (A.4|. Thus, Gronwall's Lemma yields 
the following estimate for Mj : 



Mj( X ,T) 



< 



(IrP 

\^r)\ + -^+ P (0)\M0^)\ 



+ 2W / 6XP (2R I P(y)dy\e^[\p'( S )\+T 2 d 2 ]ds; 



since $ a (Q,r) = 1, $ o (0,t) = v/T .+ d2 and < < d 2 , we have 



M^r) < |l + V2+ ^e T * [d 2 + \q\ BV ]X , (A.7) 



for |t| > <i. By a standard approximation argument we can infer that (A.7) holds for 
every q € BV(M). By ( |A6| , ( |A4| and ( |A7| , we get that 



^■(I,r) = ^(0,r)e i " i + O(i 



Again, from (A.6 1 and the above asymptotic formula, we obtain (A.l I and (A. 2 1. □ 
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In Lemmas |A.2| and |A.3| we will use the following inequality: 

. . ,, ,\ x/d 2 - 7m 
lm sm t < max < 1 , — 



2/3o 



, ie^jG [O.tt/2]. 



(A. 



Lemma A. 2. Let g be defined by (3.8 1. Then, the following asymptotic expansion 

{ x }h 



iPo({x} h -£) sini 



2(3 sin t 



p(y)dy 



O 



1 

sin t\ 



(A.9) 



holds uniformly as t — ► oo on C® for ■& 6 [0,tt/2], a; € K and £ bounded. 



Proof. The proof is a straightforward consequence of Lemma A.l and ( A. 8 1, and 
hence is omitted. □ 



Lemma A. 3. Let g be given by ( |3.8| . Then g is a bounded function of x, £ G R, 
i/£ is bounded, and t 6 $ 6 [0,7r/2]. 

Proof, (i) First, we prove an estimate for 0j(x,t 2 + cP) for |t| < ci, |Imr| < T 

and |x| < /i. By setting A = r 2 + d 2 and %(x,A) = sup |0j(s,A)| as before, from 

se[o,x] 



(A. 3 1 and since 



is bounded by the constant B = J cosh (Th) 



have 



nh 2 (Th) 
(Th) 2 ' 



we 



kl 

n s (x,X) <cosh(Tx) + B / p(j/)|x - t/|r? s (t/, X)dy, 

a 

\x\ 

Va(x,X) < B W2d\x\ + J p(y)\x-y\r] a (y,X)dy 
o 

Gronwall's Lemma yields 

d 2 h 2 



\<f>j(x,X)\ < min|cosh(T/i),\/2d/i-Blexp (b 



for Ixl < h, 1-7-1 < d and llmrl < T. 



(ii) To complete the proof, we notice that from (A. 5) it follows that 

de Th 



$j(x,t) 



< 



\r^(h,r)\ 



[V2 + dhe Th sup |0j(x,A)| 

V x£[— h,h] 



(A.10) 



and since r$j(h, t) ^ far from the poles of g, we have that ^ J |^j is bounded for 
t = /3 sini, t £ C$, $ € [0,7r/2] and for |r| < d. Thus, the assertion of the lemma 
follows from (A. 10 1 and Lemma [A. 2 □ 
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